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Abstract. Massive stars played a key role in the early evolution of the Universe. They formed 
with the first halos and started the re-ionisation. It is therefore very important to understand 
their evolution. In this review, we first recall the effect of metallicity (Z) on the evolution of 
massive stars. We then describe the strong impact of rotation induced mixing and mass loss 
at very low Z. The strong mixing leads to a significant production of primary 14 N, 13 C and 
22 Ne. Mass loss during the red supergiant stage allows the production of Wolf-Rayet stars, type 
Ib,c supernovae and possibly gamma-ray bursts (GRBs) down to almost Z = for stars more 
massive than 60 A/0. Galactic chemical evolution models calculated with models of rotating 
stars better reproduce the early evolution of N/O, C/O and 12 C/ 13 C. Finally, the impact of 
magnetic fields is discussed in the context of GRBs. 
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1. Introduction 

Massive stars started forming about 400 millions years after the Big Bang and ended 
the dark ages by re-ionising the Universe. They therefore played a key role in the early 
evolution of the Universe and it is important to understand the properties and the evolu- 
tion of the first stellar generations to determine the feedback they had on the formation 
of the first cosmic structures. It is unfortunately not possible to observe the first massive 
stars because they died a long time ago but their chemical signature can be observed in 
low mass halo stars (called EMP stars), which are so old and metal poor that the inter- 
stellar medium out of which these halo stars formed are thought to have been enriched 
by one or a few massive stars. Since the re-ionisation, massive stars have continuously 
injected kinetic energy (via various types of supernovae) and newly produced chemical 
elements (by both hydrostatic and explosive burning and s and r processes) into the 
interstellar medium of their host galaxy. They are thus important players for the chemo- 
dynamical evolution of galaxies. Most massive stars leave a remnant at their death, either 
a neutron star or a black hole, which produce pulsar or X-ray binaries. 

The evolution of stars is governed by three main parameters, which are the initial 
mass, metallicity (Z) and rotation rate. The evolution is also influenced by the presence 
of magnetic fields and of a close binary companion. For massive stars (M > 10 M©) 
around solar metallicity mass loss plays a crucial role, in some cases removing more 
than half of the initial mass. Internal mixing, induced mainly by convection and rotation 
also significantly affect the evolution of stars. In this review, after a summary of the 
properties of law-Z stars, we discuss the possible impact of rotation induced mixing and 
mass loss at low Z. We then present the implication of strong mixing and mass loss for 
the nucleosynthesis and for galactic chemical evolution in the context of extremely metal 
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poor stars. We also discuss the effects of magnetic fields. We end with conclusions and 
an outlook. 



2. Properties of non-rotating low-Z stars 



M-20 M s : metallicity effect {no rotation) 
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Figure 1. Left: H-R diagram for non-rotating 20 M Q models with Z = 10~ 8 , 10~ 5 , 0.001 and 
0.02, showing that more metal poor stars have more compact envelopes and are less likely to 
reach the red supergiant stage. Right: Central temperature versus central density diagram for 40 
Mq models with Z = 0, 10 -8 and 0.02. The evolutionary tracks start where hydrogen burning 
start and the He symbols are placed at the start of helium burning. The central conditions 
are much hotter and denser at very low Z. Note also the different initial H-burning conditions 
between Z = and 10~ 8 (explained in the text). 

The first stellar generations are different from solar metallicity stars due to their low 
metal content or absence of it. First, \aw-Z stars are more compact (see Fig. [IJ due to 
lower opacity. Second, metal free stars burn hydrogen in a core, which is denser and hotter 
due to the lack of initial CNO elements (see start of dashed curve in Fig.Q] right). This 
implies that the transition between core hydrogen and helium burning is much shorter 
and smoother. Furthermore, hydrogen burns via the pp-chain in shell burning. These 
differences make the metal free stars different from the second or later generation stars 
(|Ekstr5m et al. 2007|) . 

Third, mass loss is metallicity dependent and therefore is expected to become very 
weak at very low metallicity. The metallicity (Z) dependence of mass loss rates is usually 
described using the formula: 

M{Z) = M{Z Q ){Z/Z Q )<* (2.1) 

The exponent a varies between 0.5-0.6 (jKudritzki fc Puis 20001 [Kudritzki 2002)) and 0.7- 
0.86 (|Vink et al. 20011 IVink k de Koter 2005P for O-type and WR stars respectively (See 
Mokicm et al. 2007l for a recent comparison between mass loss prescriptions and observed 
mass loss rates). Until very recently, most models use at best the total metal content 
present at the surface of the star to determine the mass loss rate. However, the surface 
chemical composition becomes very different from the solar mixture, due either to mass 
loss in the WR stage or by internal mixing (convection and rotation) after the main 



Stellar Evolution at Low Metallicity 



103 



sequence. It is therefore important to know the contribution from each chemical species 
to opacity and mass loss. 

Recent studies (jVink et al. 20001 IVink fe de Koter 2005)) show that iron is the domi- 
nant element concerning radiation line-driven mass loss for O-type and WR stars. In the 
case of WR stars, there is however a plateau at low metallicity due to the contributions 
from light elements like carbon, nitrogen and oxygen (CNO). In between the hot and 
cool parts of the HR-diagram, mass loss is not well understood. Observations of the LBV 
stage indicate that several solar masses per year may be lost (jSmith et al. 2003)) and 
there is no indication of a metallicity dependence. In the red supergiant (RSG) stage, 
the rates generally used are still those of Nieuwenhuijzen & de Jager (1990) More recent 
observations indicate that there is a very weak dependence of dust-driven mass loss on 
metallicity and that CNO elements and especially nucleation seed components like sili- 
con and titanium are dominant (jvan Loon 20001 Ivan Loon 2006 , Fcrrar otti fe Gail~"2~0 06). 



Van Loon et al. (2005)| provide recent mass loss rate prescriptions in the RSG stage. 



In particular, the ratio of carbon to oxygen is important to determine which kind of 
molecules and dusts form. If the ratio of carbon to oxygen is larger than one, then 
carbon-rich dust would form, and more likely drive a wind since they are more opaque 
than oxygen-rich dust at low metallicity (Horner & Andersen 2007). 

Fourth, the binary interactions are probably changed by the greater compactness of 
low Z stars on the MS. Furthermore, the first generation stars below 40 Mq do not evolve 
to the RSG stage. De Mink et al. (2007) show that as Z decreases, mass transfer is more 
likely to take place after the ignition of He burning (case C) and for the very \ow-Z stars, 
which do not reach the RSG stage, only the closest binaries would still interact. 

Finally, the first stars are thought to be more massive than solar-metallicity stars 
(jBromm fc Larson 2004[ ISchneider et al. 2006)) . Note that other studies suggest that 
both very massive and low-mass stars may form at Z — ( Naka mura fc U mcmura 2001). 
Since mass loss is expected to be very low at very low metallicities, the logical deduction 
from these two arguments is that a large fraction of the first stars were very massive at 
their death (> 100 Mq) and therefore lead to the production of pair-creation supernovae 
(PCSNe). Unfortunately, the first massive stars died a long time ago and will probably 
never be detected directly (see however Scannapieco et al. 2005 ITornatore et al. 2007)) . 
There are nevertheless indirect observational constraints on the first stars coming from 
observations of the most metal-poor halo stars ([Beers fc Christlieb 2005P - These obser- 
vations do not show the peculiar chemical signature of PCSNe (strong odd-even effects 
and low zinc, see Heger & Woosley 2002 1. This probably means that at most only a 
few of these very massive stars (>100 Mq) formed or that they lost a lot of mass even 
though their initial metal content was very low as discussed in the next section. Although 
there is no signature of PCSNe at very low Z, they might occur in our local Universe 
(jSmith et a L_2007| [Larlger et al. 2007) [Woosley et al. 2007[ ). 

The topic of low Z stellar evolution is not new (see for example lChiosi 19 83 , El Ei d et al. 1983) 
ICarr et al. 198 4, Arnet t 1996)) . The observations of extremely metal-poor stars ()Beers fc C hristlieb 2005) 
have however greatly increased the interest in very metal-poor stars. There are many re- 



cent works studying the evolution of metal- free (or almost) massive (Heger fc Woosley 2002 



Limongi fc Chieffi 2005)IUmeda fc Nomoto 20 05, Mcynet et al. 2006), intermediate mass 
(|Siess et al. 2002)|Herwig 2004)ISuda et al. 2004")IGil-Pons et al. 2005P and low mass (jPicardi et al. 20041 
IWeiss et al. 2004)) stars in an attempt to explain the origin of the surface abundances ob- 



served. The fate of non-rotating massive single stars at low Z is summarised in Heger et al. (2003) 



and several groups have calculated the corresponding stellar yields (He ger fc Woosley 2002) 



Chieffi fc Limongi 2004 Tominaga et al. 2007). 
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3. Rotation, internal mixing and mass loss 

Massive star models including the effects of both mass loss and especially rotation 
better reproduce many observables around solar Z (See contributions by Meynet and 
Maeder in this volume). For example, models with rotation allow chemical surface en- 
richments already on the main sequence (MS), whereas without the inclusion of rota- 
tion, self-enrichments are only possible during the RSG stage ( |Heger fc Lan gcr 2000, 
Meynet & Maeder 2000). Rotating star models also better reproduce the WR/O ratio 
and also the ratio of type Ib+Ic to type II supernova as a function of metallicity compared 
to non-rotating models, which underestimate these ratios ( |Meynet fe Maeder 2005[ ). The 
models at very low Z presented here use the same physical ingredients as the successful 
solar Z models. The value of 300 kms -1 used as the initial rotation velocity at solar 
metallicity corresponds to an average velocity of about 220kms _1 on the main sequence 
(MS) which is close to the average observed value. See for instance |Fukuda (1982) for one 
of the first surveys and the list in Meynet 's contribution in this volume for the most recent 
surveys. It is unfortunately not possible to observe very low Z massive stars and measure 
their rotational velocity since they all died a long time ago. Higher observed ratio of Be 
to B stars in the Magellanic clouds compared to our Galaxy (jMaeder et al. 1999)) could 
point out to the fact the stars rotate faster at lower metallicities. Also a low-Z star having 
the same ratio of surface velocity to critical velocity, v (where u cr it is the velocity for 
which the centrifugal force balances the gravitational force) as a solar- Z star has a higher 
surface rotation velocity due to its smaller radius (one quarter of Zq radius for a very 
\ow-Z 20 M© star). In the models presented below, the initial ratio u/w C rit is the same or 
slightly higher than for solar Z (sec Hirs chi 20071 for more details). This corresponds to 
initial surface velocities in the range of 600 — 800 kms -1 . These fast initial rotation ve- 
locities are supported by chemical evolution models of Chiappini et al. (2006b) discussed 
in the next section. The mass loss prescriptions used in the Geneva stellar evolution code 
are described in detail in Meynet & Maeder (2005). In particular, the mass loss rates 
depend on metallicity as M ~ (Z/Zq)°- 5 , where Z is the mass fraction of heavy elements 
at the surface of the star. 

How do rotation induced processes vary with metallicity? The surface layers of massive 
stars usually accelerate due to internal transport of angular momentum from the core to 
the envelope. Since at low Z, stellar winds are weak, this angular momentum dredged 
up by meridional circulation remains in the star, and the star more easily reaches critical 
rotation. At the critical limit, matter can easily be launched into a keplerian disk which 
probably dissipates under the action of the strong radiation pressure of the star. 

The efficiency of meridional circulation (dominating the transport of angular momen- 
tum) decreases towards lower Z because the Gratton-Opik term of the vertical velocity 
of the outer cell is proportional to 1/ p. On the other hand, shear mixing (dominating the 
mixing of chemical elements) is more efficient at low Z. Indeed, the star is more compact 
and therefore the gradients of angular velocity are larger and the mixing timescale (pro- 
portional to the square of the radius) is shorter. This leads to stronger internal mixing 
of chemical elements at low Z ( |Meynet fc Maeder~ 2002). 

Figure [2] shows the evolution of the convective zones in a rotating and a non-rotating 
20 Mq models at Z = 10~ 8 . The history of convective zones (in particular the con- 
vective zones associated with shell H burning and core He burning) is strongly affected 
by rotation induced mixing. The most important rotation induced mixing takes place 
while helium is burning inside a convective core. Primary carbon and oxygen are mixed 
outside of the convective core into the H-burning shell. Once the enrichment is strong 
enough, the H-burning shell is boosted (the CNO cycle depends strongly on the carbon 
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Figure 2. Structure evolution diagram for the non-rotating (left) and rotating (right) 20 Mq 
models at Z — 10 -8 . Coloured areas correspond to convective zones along the Lagrangian mass 
coordinate as a function of the time left until the core collapse. The burning stage abbreviations 
are given below the time axis. Rotation strongly affects shell H burning and core He burning. 



and oxygen mixing at such low initial metallicities). The shell then becomes convective 
and leads to an important primary nitrogen production. In response to the shell boost, 
the core expands and the convective core mass decreases. At the end of He burning, the 
CO core is less massive than in the non-rotating model (see Fig. \2§ . Additional convective 
and rotational mixing brings the primary CNO to the surface of the star. This has con- 
sequences for the stellar yields. The yield of 16 being closely correlated with the mass 
of the CO core, it is therefore reduced due to the strong mixing. At the same time the 
carbon yield is slightly increased. The relatively "low" oxygen yields and "high" carbon 
yields are produced over a large mass range at Z = 10~ 8 (Hirschi 2007). This could be 
an explanation for the possible high [C/O] ratio observed in the most metal-poor halo 
stars (ratio between the surface abundances of carbon and oxygen relative to solar; see 
Fig. 14 in |Spite et al. 20050 . 

Models of metal- free stars including the effect of rotation (Ekstrom et al. 2006) show 
that stars may lose up to 10 % of their initial mass due to the star rotating at its 
critical limit (also called break-up limit). The mass loss due to the star reaching the 
critical limit is non-negligible but not important enough to change drastically the fate of 
the metal-free stars. The situation is very different at very low but non-zero metallicity 
( |Meynet et al. 2006| IHirschi 2007| . The total mass of an 85 Mq model at Z = 10 is 
shown in Fig. [3] by the top solid line. This model, like metal-free models, loses around 
5% of its initial mass when its surface reaches break-up velocities in the second part of 
the MS. At the end of core H burning, the core contracts and the envelope expands, 
thus decreasing the surface velocity and its ratio to the critical velocity. The mass loss 
rate becomes very low again until the star crosses the HR diagram and reaches the RSG 
stage. In the cooler part of the H-R diagram, the mass loss becomes very important. 
This is due to the dredge-up by the convective envelope of CNO elements to the surface 
increasing its overall metallicity. The total metallicity, Z, is used in this model (including 
CNO elements) for the metallicity dependence of the mass loss. Therefore depending on 
how much CNO is brought up to the surface, the mass loss becomes very large again. The 
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CNO brought to the surface comes from primary C and O produced in the He-burning 
region and from primary N produced in the H-burning one. 

M= 85 M e , v nl = 800 km s" 1 , Z = le-08 
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Figure 3. Structure evolution diagram (same as Fig. 2) for a 85 Mq model with u; n j = 800 
kms -1 at Z — 10 -8 . The top solid line shows the total mass of the star. A strong mass loss 
during the RSG stage removes a large fraction of total mass of the star. 

Could such low- .2? stars undergo dust-driven winds? For this to occur, the surface ef- 
fective temperature needs to be low enough (usually log(T e g) < 3.6) and carbon needs 
to be more abundant than oxygen. This last condition is fulfilled in our 85 Mq model. 
However, it is presently unclear if: 1) Extremely low-Z stars reach such low effective tem- 
peratures. This depends on the opacity and the opacity tables used in our calculations did 
not account for the non-standard mixture of metals (high CNO and low iron abundance, 
sec Marigo ~2002| for possible effects). 2) At such low Z, enough metal is present to allow 
dust formation. Indeed, nucleation seeds (probably involving titanium) are necessary to 
form C-rich dust. There may also be other important types of wind, like chromospheric 
activity-driven, pulsation-driven, thermally-driven or continuum-driven winds. 

The fate of rotating stars at very low Z is therefore probably the following: 

• M < 40 M Q : Mass loss is insignificant and matter is only ejected into the ISM 
during the SN explosion (see contributions by Nomoto et al, Limongi et al and Frohlich 
et al in this volume), which could be very energetic if fast rotation is still present in the 
core at the core collapse. 

• 40 Mq < M < 60 Mq : Mass loss (at critical rotation and in the RSG stage) removes 
10-20% of the initial mass of the star. The star probably dies as a black hole without a 
SN explosion and therefore the feedback into the ISM is only due to stellar winds, which 
are slow. 

• M > 60 Mq : A strong mass loss removes a significant amount of mass and the stars 
enter the WR phase. These stars therefore die as type Ib/c SNe and possibly as GRBs. 



4. Nucleosynthesis and galactic chemical evolution 

Rotation induced mixing leads to the production of primary nitrogen, 13 C and 22 Ne. In 
this section, we compare the chemical composition of our models with carbon-rich EMP 
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stars and include our stellar yields in a galactic chemical evolution (GCE) model and 
compare the GCE model with observations of EMP stars. 

4.1. The most metal-poor star known to date, HE1327-2326 

Wind composition, [Fe/H] — 6.6 (Z=10- 8 ) 
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Figure 4. Composition in [X/Fe] of the stellar wind for the Z — 10 -8 models (solid lines). 
For HE1327-2326 (red stars), the best fit for the CNO elements is obtained by diluting the 
composition of the wind of the 40 Mq model by a factor 600 (see Hirschi 2007 for more details) . 




Significant mass loss in very low-Z massive stars offers an interesting explanation for 
the strong enrichment in CNO elements of the most metal-poor stars observed in the halo 
of the galaxy (see |Meynet et al. 2006 [Hirschi 2007[) . The most metal-poor star known 
to date, HE1327-2326 (|Frebel et al. 2006[) is characterised by very high N, C and O 
abundances, high Na, Mg and Al abundances, a weak s-process enrichment and depleted 
lithium. The star is not evolved so has not had time to bring self-produced CNO elements 
to its surface and is most likely a subgiant. By using one or a few SNe and using a very 
large mass cut, Limongi et al. (2003) and Iwamoto et al. (2005) are able to reproduce 
the abundance of most elements. However they are not able to reproduce the nitrogen 
surface abundance of HE1327-2326 without rotational mixing. The abundance pattern 
observed at the surface of that star present many similarities with the abundance pattern 
obtained in the winds of very metal poor fast rotating massive star models. HE1327- 
2326 may therefore have formed from gas, which was mainly enriched by stellar winds of 
rotating very low metallicity stars. In this scenario, a first generation of stars (PopIII) 
pollutes the interstellar medium to very low metallicities ([Fe/H]~-6). Then a PopII.5 
star (Hirschi 2005) like the 40 M© model calculated here pollutes (mainly through its 
wind) the interstellar medium out of which HE1327-2326 forms. This would mean that 
HE1327-2326 is a third generation star. In this scenario, the CNO abundances are well 
reproduced, in particular that of nitrogen, which according to the latest values for a 
subgiant (see lFrebel et al. 2006[) is 0.9 dex higher in [X/Fe] than oxygen. This is shown 
in Fig. 2] where the abundances of HE1327-2326 are represented by the red stars and 
the best fit is obtained by diluting the composition of the wind of the 40 M© model by 
a factor 600. When the SN contribution is added, the [X/Fe] ratio is usually lower for 
nitrogen than for oxygen. It is interesting to note that the very high CNO yields of the 
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Figure 5. Chemical evolution model predictions of the N/O and C/O evolution, in the galactic 
halo, for different stellar evolution inputs. The solid curves show the predictions of a model with- 
out fast rotators at low metallicities. The dashed lines show the effect of includin g a population 
of fast rotators at low metallicities. For the data see Chiappini et al. (2006b) and references 
therein. 

40 M & stars brings the total metallicity Z above the limit for low mass star formation 



obtained in Bromm & Loeb (2003) 



4.2. Primary nitrogen and 13 C 

The high N/O plateau values observed at the surface of very metal poor halo stars 
require very efficient sources of primary nitrogen. Rotating massive stars can inject in 
a short timsecale large amount of primary N at low Z. They are therefore very good 
candidates to explain the N/O plateau observed at very low metallicity. According to 
the heuristic model of Chiappini et al. (2005) a primary nitrogen production of about 
0.15 Mq per star is necessary. Upon the inclusion of the stellar yields including the 
effects of fast rotation at Z=10 -8 in a chemical evolution model for the galactic halo with 
infall and outflow, both high N/O and C/O ratios are obtained in the very metal-poor 
metallicity range in agreement with observations (see details in Chiappi ni et al. 2006a[ ). 
This model is shown in Fig. [5] (dashed magenta curve). In the same figure, a model 
computed without fast rotators (solid black curve) is also shown. Fast rotation enhances 
the nitrogen production by ^3 orders of magnitude. These results also offer a natural 
explanation for the large scatter observed in the N/O abundance ratio of normal metal- 
poor halo stars: given the strong dependency of the nitrogen yields on the rotational 
velocity of the star, we expect a scatter in the N/O ratio which could be the consequence 
of the distribution of the stellar rotational velocities as a function of metallicity. 

As explained above, the strong production of primary nitrogen is linked to a very active 
H-burning shell and therefore a smaller helium core. As a consequence, less carbon is 
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turned into oxygen, producing high C/O ratios. Although the abundance data for C/O 
is still very uncertain, a C/O upturn at low metallicities is suggested by observations 
(see |Asplund 2005| and references therein). Note that this upturn is now also observed in 
very metal poor DLA systems (see the paper by M. Pettini in this volume). 

In addition, stellar models of fast rotators have a great impact on the evolution of the 
12 C/ 13 C ratio at very low metallicities (Chiappini et al. 2007). In this case, we predict 
that, if fast rotating massive stars were common phenomena in the early Universe, the 
primordial interstellar medium of galaxies with a star formation history similar to the 
one inferred for our galactic halo should have 12 C/ 13 C ratios between 30-300. Without 
fast rotators, the predicted 12 C/ 13 C ratios would be ~ 4500 at [Fe/H] = —3.5, increasing 
to - 31000 at around [Fe/H] = -5.0 (see Fig. 0. 

Current data on EMP giant normal stars in the galactic halo ( Spit e et al. 2006| ) agree 
better with chemical evolution models including fast rotators. The expected difference in 
the 12 C/ 13 C ratios, after accounting for the effects of the first dredge-up (indicated by the 
arrows in Fig. [6|), between our predictions with/without fast rotators is of the order of a 
factor of 2-3. However, larger differences (a factor of ~ 60 — 90) are expected for giants at 
[Fe/H]= —5 or turnoff stars already at [Fe/H]= —3.5. To test our predictions, challenging 
measurements of the 12 C/ 13 C in more extremely metal-poor giants and turnoff stars are 
required. 



4.3. Primary 22 Ne and s process at low Z 

Models at Z = 10~ 8 show a production of primary 22 Ne during He burning. We also 
started calculating models at different Z to determine over which Z range the primary 
production of 22 Ne and also 14 N is important. In Fig. \7\ we show the properties of a 20 
M© model at Z — 10~ 6 up to the end of He burning. Around 0.5% (in mass fraction) of 
22 Ne is burnt during core He burning and therefore leads to a significant neutron release. 
Studies are underway to determine how much s process can be produced in these models. 
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Figure 7. Left: Structure evolution diagram (see description of Fig. [2} for a rotating 20 Mq 
models at Z — 10 -6 during H- and He-burning phases. Right: Chemical composition at the 
end of core He burning. Just above the core, one sees that the maximum abundance of 22 Ne is 
around 1% in mass fraction and at the end of core He burning, around 0.5% is burnt in the core, 
providing plenty of neutrons for s process. 



5. Magnetic fields and GRBs 

In this last section, we discuss the impact of magnetic fields. Models of rotating stars, 
which do not include the effect of magnetic fields predict gamma-ray bursts (GRBs) at 
almost all Z (Hirschi et al. 2005, Hirschi 2007). However, they also overestimate the ini- 
tial rotation rate of pulsars. The inclusion of the effects of magnetic fields according to 



Spruit (2002) allows a better reproduction of the initial pulsar periods ( |Heger et al. 2005 ). 



Along with gravity waves, magnetic fields are also one possible cause for the fiat rotation 
profile of the Sun ([Eggcnbcrger et al. 2005). Although it becomes much harder for the 
core to retain enough angular momentum until the core collapse, there is still an evolu- 
tionary scenario, the so-called chemically homogeneous evolution, leading to the produc- 
tion of fast rotating cores at the pre-SN stage and therefore enabling MHD explosions 
and GRBs (see Yoon et al. 2006, |Woosley fc Heger 2006| and also Yoon's contribution 
in this volume). The theoretical GRB event rates obtained by Yoon et al. (2006) are in 



good agreement with observations apart from the upper metallicity limit, which is lower 



than the observed one (see contribution by Stanek in this volume). Yoon et al. (2006) 
also predict that at Z = 10 -5 , a large fraction of massive stars arc GRB progenitors. 
We have calculated 40 M Q models at Z — 10~ 5 with Wmi/w cr it = 0.59 and at Z = 10~ 8 
with Vini/t>crit = 0.55. The model at Z = 10~ 5 confirms the possibility of producing 
GRBs down to very low Z. However, the model at Z = 10~ 8 does not rotate fast enough 
to evolve chemically homogeneously. The difficulty of the very low Z models to evolve 
chemically homogeneously is due to the weakening of the meridional circulation. Indeed, 
in models including magnetic fields, meridional circulation becomes the dominant term 
for the mixing of chemical species (see Maeder & Meynet 2005). The meridional circu- 
lation becomes weaker at low Z because the meridional currents are less efficient in a 
denser medium, which is the case since low-Z stars are more compact. This means that 
not all stars in the first stellar generations will produce GRBs in this way. Finally, it is 
interesting to note that the presence of magnetic fields in metal-free stars may enhance 
mass loss significantly (see contribution by Ekstrom et al in this volume). 
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6. Conclusions and outlook 

The inclusion of the effects of rotation changes significantly the simple picture in 
which stellar evolution at low Z is just stellar evolution without mass loss. A strong 
mixing is induced between the helium and hydrogen burning layers leading to a significant 
production of primary 14 N, 13 C and 22 Ne. Rotating stellar models also predict a strong 
mass loss during the RSG stage for stars more massive than 60 M . The chemical 
composition of the stellar winds is compatible with the CNO abundance observed in the 
most metal-poor star known to date, HE1327-2326. GCE models including the stellar 
yields of these rotating star models are able to better reproduce the early evolution of 
N/O, C/O and 12 C/ 13 C in our galaxy. These models predict a large neutron release during 
core He burning and thus a strong possibility of an s process at very low Z. These models 
predict the formation of WR and type Ib/c SNe down to almost Z=0. The inclusion of 
magnetic fields slows down the core of the stars and therefore reduces the probability 
of producing GRBs at metallicities around that of the Magellanic Clouds but GRBs are 
still predicted from single star models down to very low Z. 

Large surveys of EMP stars (SEGUE, LAMOST, and OZ surveys), of GRBs and SNe 
(Swift and GLAST satellites) and of massive stars (e. g. VLT FLAMES survey) are 
underway and will bring more information and constraints on the evolution of massive 
stars at low Z. On the theoretical side, more models are necessary to fully understand 
and study the complex interplay between rotation, magnetic fields, mass loss and binary 
interactions at different metallicities. Large grids of models at low Z will have many 
applications, for example to study the evolution of massive stars and their feedback in 
high redshift objects like Lyman-break galaxies and damped Ly-alpha systems. 
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